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DETERMINING THE COMPOSITION OF PHYSICAL 
STRUCTURES FROM BOUNDARY MEASUREMENTS, 
INCLUDING THE COMPOSITION OF SUBSCRIBER LOOPS 

CROSS-REFERENCE TO RELATED APPLICATION 

This application is a non-provisional application of provisional application 
Serial No. 60/158,156 filed October 8, 1999. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

This invention relates to physical systems identifiable via configuration 
parameters and, more particularly, to methodologies and concomitant systems for 
determining the composition of the physical system via measurements at system 
boundaries, including determining the composition of subscriber loops. 

2. Description of the Background Art 

There has been considerable recent interest in the development of 
techniques to determine the composition of physical structures from measurements of the 
physical structure at boundary points of the structure. Representative of the technology 
applied to the imaging of objects, such as tumors or masses, in the human body, is the 
technique reported inU. S. Patent No. 5,747,810 issued to Schotland and entitled 
"Simultaneous Absorption and Diffusion Tomography System and Method Using Direct 
Reconstruction of Scattered Radiation." As disclosed, the method for the direct 
reconstruction of an object from measurements of the transmitted intensity of diffusively 
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scattered radiation is effected by irradiating the object with a source of continuous wave 
radiation. The transmitted intensity is related to the absorption and the diffusion 
coefficients of the object, which effectively determine or define the object, by an integral 
operator. The image of the object is directly reconstructed by executing a prescribed 
mathematical algorithm, as determined with reference to the integral operator, on the 
transmitted intensity of the diffusively scattered radiation. 

With respect another class of physical structures, namely telephone 
subscriber loops, wherein a loop typically connects a customer with a local telephone 
central office and is composed of lengths of cable such as, for example 26 gauge or 24 
gauge cable, there has been considerable interest shown in devising a technique to 
determine the composition of a loop from so-called single-ended measurements in order to 
qualify such loops for high-speed digital transmission. It is especially desirable to estimate 
the configuration of a subscriber loop from measurements made at the input of the loop at 
the central office. For example, one might measure the complex input impedance of the 
loop in the frequency domain or the time-domain echo at the input to the loop. From 
these measurements, the composition of the loop is estimated using identifiable 
characteristics in the response, such as peaks of return signals or time intervals between 
peaks in the response signal (generally referred to as "time domain reflectometer" 
approaches). Moreover, based upon estimates of the loop configuration from these 
measurements, it is further possible to estimate the transmission characteristics of the loop 
to the customer end. 
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SUMMARY OF THE INVENTION 

These shortcomings, as well as other limitations and deficiencies are 
obviated, in accordance with the present invention relating to the composition of physical 
structures, for which a subscriber loop is representative, by devising explicit inversion 
formulae for determining the composition of the physical structure in terms of a parameter 
or configuration space. This approach, when applied to loops, allows the separation of 
loop make-up into parameters such as length of a cable section and its corresponding 
gauge. 

In accordance with a broad method aspect of the present invention, a 
method for determining configuration parameters describing a physical system includes: 
(a) measuring an output signal from the system in response to an input signal, the output 
signal being related to the configuration parameters by a linear operator; and (b) directly 
reconstructing each of the configuration parameters by applying a prescribed mathematical 
algorithm to the output signal. 

The broad system aspect of the present invention is commensurate with this 
broad method aspect. 
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BRIEF DESCRIPITON OF THE DRAWING 

The teachings of the present invention can be readily understood by 
considering the following detailed description in conjunction with the accompanying 
drawings, in which: 

FIG. 1 illustrates a high-level block diagram of system equipment to 
directly reconstruct configuration parameters, 

FIG. 2 illustrates an operational flow diagram to generate data for the 
methodology of the present invention, and 

FIG, 3 illustrates an operational flow diagram to generate data for the 
methodology of the present invention as applied to subscriber loops. 

DETAILED DESCRIPTION 
Function Theoretic Basis 

Let x be a vector that defines the loop configuration; its components are 
the lengths of the various nominally uniform segments of the loop, such as main sections 
of constant gauge and no bridged tap, bridged taps of constant gauge, and so forth. There 
may be other types of quantities that are dimensions in the configuration space as well, 
such as the real part of a terminating impedance at a specified frequency. Some of the 
components of the configuration vector may be quantized, e.g., wire gauge. The complex 

input impedance of the loop might then be represented by x) . The following is a 
method of making estimates of x from /zf <2>,x) , whatever it represents, in the form 
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x = \h(a) y x)g{6))d6) (1) 

Q 

where x is the estimated configuration vector, g{a)) is a function that is determined herein 
and Q is the domain over which measurements are made. 

The technique that is illustrated for frequency domain measurements is then 
generalized to a general class of near-end measurements that includes time-domain 
measurements. Approximation methods are also discussed. In addition, it is recognized 
that the technique revealed is generally applicable to a wide variety of system identification 
problems. 

2. ESTIMATING THE CONFIGURATION USING MEASUREMENTS IN THE 
FREQUENCY DOMAIN. 

We consider the response function, x) , as the kernel of an integral 

operator that maps functions on the space of configuration vectors onto functions of the 
frequency. Thus, we can write a transform of the form 

A 

f(co)=lh*(a>,x)f{x)dX (2) 

X 

where the asterisk indicates the complex conjugate, dXis an elemental volume in 

A 

configuration space, \dX is an integral over the whole configuration space, and / (#>)is 

x 

the transform of f(x) . Now, consider the integral operator with kernel K{y,co) defined 
by 

K{v y m)=\h{v^)h\(D^)dX (3) 
x 
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This is a positive definite kernel because there can be no function /(x) that has a zero 
transform because of causality. Consequently, the kernel can be inverted on the class of 
functions given generically by equation (2), and the result may be applied to equation (2) 
to yield 

A 

/(v)= lKT\v,m)]l{m$f{^da (4) 

Q X 

with K~ l (v,a>) is the kernel of the integral operator that is the inverse of the operator 
with kernel K(v,(d), and there are some other some obvious definitions. The domain of 
integration in equation (4), viz., Q, might, for example, be a convenient pass-band of 
frequencies, symmetric about zero but with a gap that includes 0 (DC) on the real 
frequency line. A pass-band might be chosen to avoid the low-frequency distortion that is 
inherent in copper cables. The integral operator with kernel K(v,to)\s compact, so that 
its spectrum has a limit point at zero. Consequently, the inverse may have to be 
generalized to avoid the exigencies of measurement imprecision and noise. 

If we now multiply both sides of equation (4) by h(v,y) we get 

f(y) = J J jA*(v,y)^-Kv^W^^)/(^)^^ (5) 

QQX 

because the operator that acts on f(x) on the right is a projection operator. 
To apply equation (5) to the case at hand suppose, for example, that f{x) is the 
configuration vector, x , and that h{v,y) is the result of measurement on a loop with 
configuration vector y . Then equation (5) yields 
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y=\h*(v,y)g(v)dn (6) 
n 

where y is an estimate of the configuration vector and 

g(v) = \ \k-'(v, co)h(a),x)xdXda (7) 

xa 

Equation (6) is interpreted as follows: the configuration vector can be found by taking the 
inner product of the results of the frequency domain measurements with the predetermined 
function g(v) . A discussion of the relationship between y and y is given in the sequel. 

Application of this result is illustrated FIG. 1. In the system of FIG. 1, 
subscriber loop 101 is driven by swept frequency source 110. (Source 1 10 may actually 
be a sequence of measurements that are made at a discrete set of frequencies, in which 
case, the integrations in the above that are over frequency should be replaced by 
summations.) Source 110 might, for example, be a controlled current source; the 

response or output signal, shown as /*(v,y) , might be the input voltage to the line. In 
effect, the measurement setup determines the complex driving point impedance at the 
input to the loop. 

The functions gjc(v) would either be determined in advance of 
measurement, or could be calculated and stored in a processor (not shown) at the time of 
measurement using stored characteristic transmission data on wire-pair cables. 
Alternatively, the test system might be linked to a back-end or remote computer that 
would conduct such calculations and then relay results to the test system. In either event, 
interpolation might be needed between stored or computed values of gjc(v) to obtain 
values that can be used with the measured values of the system response. 
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Integration is carried out by processors 121, 122, . . 123 to provide the 
estimates of the configuration parameters of vector y . 

3. GENERALIZATION OF THE ESTIMATION TECHNIQUE 

The technique that was elucidated in section (2) can, in an obvious way, 
also be applied to measurements of responses that are made in the time domain. All that is 
required is to think of the variables considered as frequencies in section (2) should be 
considered as times and to make other obvious adjustments, mutatis mutandis. In fact, the 
results of section (2) can be made generic to any type of measurement that will represent a 
one-to-one map of the configuration space. 

To the end of developing a generic formula for estimation of a loop 
configuration vector, define h(x) as the vector or function quantity that is measured when 
the loop configuration vector is x . Then, define the operator Hon all scalar functions on 
the space of configuration vectors, X y by 

Hf = Jh(x)/(x)d3T (8) 

X 

where f is a vector notation for the scalar function on X . Note that f may be 
generalized in an obvious way to have several scalar-function components. To exemplify, 
the configuration vector itself is one such multi-component function. 
Now, form the operator 

P = H(Hfi) *H (9) 

As in section (2), a generalized inverse of HH may be required. P is obviously a 
projection operator, (because it is idempotent and symmetric) so that for any f e f(X) , 

-8- 
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Pf = f (10) 
In equation (10), we tacitly assume that f lies in the range of the projection operator. 
Consequently, if we wanted to use the measurements of the loop, i.e., h(x) to determine 
the value of an arbitrary Sanction, say the scalar function /(X) , for the loop, we 
would a priori calculate 

g/ =(HH) _1 Hf (11) 



where, in this instance, f = f(X) . Then we would use the measurement of the loop, viz., 
h(x) , to calculate /(x) , using equation (9): 

/(x) = (h(x), g/ ) = H(HH) _1 Hf (12) 

What equation (12) does for a scalar function can, as indicated above, be done for several 
scalar functions so that vector functions can be resolved in this way. In particular, the 
vector function x can be treated in this way since it is a finite-dimensional function in the 
configuration space. Thus, we can define a vector function g x by 

g x =(HH) _1 Hx (13) 

Then, if the measurement result is h(y) , as it would be on a loop with configuration 
vector y , the estimate of the configuration vector would be 

* = (h(y),g x ) = H(y)(Hfi) X Hx (14) 

It should be noted that the technique that has been described above and is described in the 
sequel is applicable to linear operators Hthat are more general than the integral operator 
given in equation (8) and, in fact, the technique described herein is more general than 
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resolution of a telephone subscriber's loop configuration. Thus, what is described is 
applicable to a wide variety of situations where measurements are made on a system that 
depends on the system configuration and it is desired to estimate the configuration of the 
system from the measurements. 

4. DETERMINING THE INVERSE APPEARING IN (14) 

One way to determine the inverse of HH uses a singular value 
decomposition, i.e., solve the eigenvalue problem 



Then, when operating on vector quantities of the form Hf , HH is non-singular and has 
the generic form 



The summation in equation (16) may be truncated, as appropriate for practical 



We can also use the singular value decomposition to represent H and its transpose. Thus, 



(15) 




(16) 



computation purposes. Also, this is not the only way that HH can be inverted. 



if 



(17) 



it is apparent that 




(18) 



so that we can write 
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Again, truncating the summation may be used as an approximation, in accord with using a 
generalized inverse of HH . 

Another way of determining the inverse of the operator might be simply to 
represent the operator in some finite dimensional function space as a matrix. Then, the 
inverse could be found using a computer routine that inverts the matrix. Using 
measurements made at a discrete set of frequencies, as would be done in any computer 
program, would be an example of such an approach. 

5. SOME APPLICATIONS 

A. TIME DOMAIN MEASUREMENTS 

One could use the results of a time-domain reflectometer measurement in the above 
procedure. Another way to use time domain measurements is by means of an adaptive 
echo canceller at the input to the loop. The echo canceller settings that the canceller 
adapted to would be the data that is used in the above. For the purposes of illustration, 
we will assume here that the time-continuous echo is available as a response, so that the 
response function is 

h(x) = «(/,*) (20) 

where e(t,x) is the echo, at time t , from a loop with configuration vector x , when the 
loop is driven by a pulse of prescribed shape, e.g., a raised cosine pulse. The space that 

-li - 
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the response lies in is now the space of square integrable functions on the positive, semi- 
infinite, time line. The inner product between two elements, f andg, denoted by (f ,g) , is 
given by 

(f,g) = [f(t)g(t)dt (21) 

Note that the upper limit on the integration in (21) is infinite; in practice, this would be a 
finite upper limit, making the dimension of the space of functions considered denumerably 
infinite. Furthermore, from a computational point of view the function space must be 
finite-dimensional. 

The operator His given by 

H/(x)=M*,x)/(x)^ (22) 

x 

so that Hf may be viewed as a superposition of infinitesimal echoes {/(Jf)^} . The 

operator HH is just an integral operator with domain and range the space of functions that 
are square-integrable on the positive semi-infinite time line, and with kernel 

K(t,T) = \e(t,MT>x)dX (23) 

x 

If we denote the kernel of the generalized inverse of this operator by K~ l (t, r) , equation 
(23) becomes 

00 

g* = &(') = J \K l {t,r)e{r,^dXdr (24) 

0 X 

-12- 
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and the estimated value of the configuration vector, given measurement e(t,y) , would be 

00 

M«(*,y)&<f)<* (25) 

0 

It should be noted here that e(*,x) would be a scalar function in most applications, but 
that g x (t) is a finite-dimensional vector function of time, so that the configuration estimate 
is of the same dimension as e(t,x) . 

B. Sampled Configuration Space 

In most applications it is not necessary to know the system configuration exactly so that 
the configuration space can be sampled or quantized. This does not change the 
fundamental character of the above considerations. The configuration sample values 
generate a finite collection of response functions that span the response space so that the 
general formulation of section (3) applies. The necessary modifications are only that 
integral operator become matrices and integrations become summations with a finite 
number of terms. However, it may be necessary to use generalized inverses of the 
involved matrices, possibly reducing the effective dimensionality of the measurement space 
and limiting the dimensionality of the possible configuration space. 

C. Difference Responses 

It may frequently happen that an approximation to the loop configuration 
can be made a priori. For example, a preliminary measurement may indicate that the loop 
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is very long or that there is a bridged tap at some point on the line that is masking, to 
some extent, the effects of the configuration beyond the bridged tap. In this case, the 
response used in the above considerations can be the difference between a measurement of 
the actual loop and the reference measurement. This is, of course, similar to the idea of 
using a bridge to make measurements. 

In another situation, we might expect that the loop configuration lies within 
a certain subset, X Q , of the overall configuration space. E. g., an initial loop measurement 
using a technique other than that revealed herein, might give an approximate volume in 
configuration space within which the loop configuration vector lies. We might then use 
that subset for the configuration space, rather than the whole configuration space, in the 
above discussed calculations, with a reference response being that of a loop with a 
configuration somewhere in the middle of the subset. The response function then used 
would be the difference between the actual measurement and the measurement that the 
reference configuration would yield. 

Hence, suppose that there is a reference response, designated by h 0 and 
that the actual result of loop measurement is h x (x) . Then, the response that would be 
used in the above calculations would be 

h(x) = h 1 W-h. (26) 

The theory developed above could then be applied using the h(x) on the left in equation 
(26). In fact, the operator H can now be written in the form 

Hf = J hito/WdX- - h 0 \f(*)dX (27) 

x 0 x 0 



-14- 
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so that the h 0 part of the H operator has the effect of subtracting out the effect of the 
average of any function over the domain of integration. It must be remembered that the 
domain of integration in equation (27) might not be the whole set of possible 
configurations, but a smaller set, containing h 0 , as ascertained by preliminary 
considerations. 

Now, the operator HH takes the form 

HH = J [h^x) - hjhjd) - h 0 }lX (28) 



If we choose h 0 so that it is the average of hj(x) over the configuration set: 



hJdX=\h l (x)dX (29) 

X 0 x. 



Equation (28) assumes the simple form: 



HH = J h&hjixtfX - h 0 h 0 \dX (30) 

x 0 x 0 



In addition, the choice in equation (29) converts equation (27) to the form 



Hf = fh 1 (x)/(x) a HT-h /° = Jh^/WflK-- h 0 (/) (31) 

X. J OA X, 

X 0 

where (/ ) is the average of the function /(x) over the set of configurations deemed 
feasible (X 0 ). 

-15- 
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In the measurement process, a configuration specified by a vector y 
would, of course, have a result h t (y) , and h G would have to be subtracted from the 
measurements in processing the data. 

E. Iteration 

It should be noted that this idea can be used as the basis of an iteration 
process that successively refines the estimate at each stage of the iteration. For example, 
one might use a smaller subset of the overall configuration space at each step of the 
iteration, with h 0 chosen in the middle of the progressively smaller subsets. The iteration 
process could be terminated when the norm of the difference between the measurement 
and the reference response was smaller than some predetermined value. Of course, this 
does not preclude the use of other iteration methods such as steepest descent methods. 

F. Projection Annihilation 

It may sometimes happen that the range of the projection operator given in 
equation (9) is a proper subspace of the domain of the projection. This might happen, for 
example, if generalized inverses of operators are used, as alluded to above. In this case 
some functions on the configuration space will be annihilated by the projection operator. 
Consequently, projecting the configuration vector function will result in a function on the 
configuration space that is not equal to the configuration vector, but a function of the 
configuration vector. This does not alter the feasibility of the above procedure as long as 
the dimensionality of the measurement space, i.e., the dimensionality of the vector 
measured, is greater than or equal to the dimensionality of the configuration space (not the 

-16- 
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dimensionality of the space of square-integrable functions on the configuration space). 
Suppose, for example that the dimensionality of the measurement space is equal to the 
dimensionality of the putative configuration space. Then, by the inverse function theorem 
the measurements can be inverted to determine the configuration vector under fairly 
general conditions. All that is needed is, in effect, a generalized graph of the 
measurements as a function of the configuration vector, which would serve as a look-up 
table at the time of measurement. To determine such a look up table it is only necessary 
calculate h(y) for each configuration vector y prior to measurement. Of course, only a 
finite number of such responses could be stored. However, interpolation could be used at 
the time of measurement. 

G. Measurement Noise 

Some noise is present in all measurements. The above procedures can 
easily be modified so that the result of processing is a best estimate of any function of the 
configuration vector in the presence of noise. This technique will now be outlined. 
Suppose that there is noise at the time of measurement so that what is measured is not 
h(x) but h(x) + n where n is a random noise vector with covariance operator N, defined 
by 

((x,nXy 5 n))-(x ? Ny) (31) 
for all x and y in response space. In equation (3 1), < > indicates expectation in the noise 
probability space and ( ) indicates inner product in response space. Next, let gf be a 
response space vector that estimates the value of a function, /(x) , on the configuration 

- 17 - 
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space by linear filtering operation, so that the estimate of f(x) , which we will denote by 
is given by 

/k) = ( g/ ,h( X ) + ii) (32) 
The mean-square error in this estimate, integrated over configuration space, is 

« = \z<^fM + n)- /(x)][( g/ ,h(x) + n) -/(x)] > dX (33) 

or, 

e = (g / ,Kg / ) + (g p Ng / )j^-2fg / Jh( X )/( X )cff) + J/ 2 ( X )^ (34) 

X ^ X * X 

Varying % f to minimize this error leads to 

J(h(x)h(x) + N)^T g/ = lh(x)f(x)dX (35) 

X X 

or, with an obvious definition, 

(K + N) \dXg f = f h(x)f(x)dX (36) 

x x 

Thus, the introduction of noise into the measurement process modifies the procedures 
given above by adding the noise covariance matrix to the operator previously used to 
obtain g f . This has an advantage when the K operator is compact so that its eigenvalues 

have zero as a limit point, as explained above. If the noise is white, the spectrum of the 
noise covariance operator is flat, so that the limit point of the spectrum of the operator 
K + N is strictly greater than zero, thus allowing the operator K + N to be inverted on the 
whole of the response space. 

-18- 
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This suggests a method of iteration in the configuration resolution process. 
One can begin to resolve the configuration by assuming a relatively large amount of noise 
(even if no appreciable noise is actually present) so that relatively imprecise resolution is 
obtainable. This is tantamount to assuming a coarse-grained resolution, permitting rough 
estimation of the configuration vector as a first step. The putative noise might, in fact, be 
considered as the quantization noise inherent in the granularity of the resolution. Then, 
after the initial estimates using the technique of this section, the subset of configuration 
space that is considered possible can be restricted. Then, the putative noise level can be 
reduced to reflect the smaller domain of consideration and the method of this section can 
be reapplied. Continuation of this process leads to more and more refined estimates of the 
system configuration. 

Operational Flow 

The flow diagram of FIG. 2 depicts the steps of the operational flow presuming frequency 
measurements are made, as follows: 

FLOW SEQUENCE FOR FREQUENCY MEASUREMENTS 

Step 205 — Define the subset of loop configuration space that is of interest, that is, 
determine which loop configurations are expected to be possible. This includes selecting 
which dimensions in the configuration space are relevant (cable lengths, bridged tap 
lengths, gauges, etc.) 
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Step 210 — Select a finite collection of configurations within this subset so that the 
response for any configuration within the subset can be adequately approximated by a 
response from the selected finite collection. 

Step 215 — Determine the frequency responses, e.g., magnitude of driving point 
impedances at the loop input, of the configurations within the finite collection. 

Step 220 - Store the arrays representing the frequency responses. 

Step 225 — Form an array of all of the pairs of frequencies that can be formed from the 
set of frequencies for which measurements will be made. 

Step 230 — For each pair of frequencies stored in step 225, multiply the frequency 
responses for each configuration in the collection from step 210. Sum these products over 
all configurations from step 210. Store the result as a matrix. 

Step 235 — Invert the matrix found in step 230. The elements of this inverted matrix 
are indexed by pairs of frequencies, which are referred to as the input and output 
frequencies. 

Step 240 — For each frequency response determined in step 215, multiply the 
frequency response by the configuration vector that the response represents. Sum these 
products over all configuration vectors representing the configurations selected in step 
210. 

Step 245 — Multiply each element of the matrix found in step 235 by the result of step 
240 with the same frequency as the output frequency of the matrix element. 

Step 250 Sum the result of step 235 over all output frequencies and for each input 
frequency. 
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Step 255 Store the result of step 250. Call the stored vector function of frequency 
the g-function array. 

Step 260 ~ Measure the loop to determine the measured frequency response. 

Step 265 « For each frequency at which measurements are made, multiply the result of 
measurement by the vector in the g-function array. 

Step 270 — Sum the results of step 265 over all measurement frequencies. 

Operational Flow Example in Terms of General Subscriber Loop Parameters 

The flow diagram in FIG. 3 depicts the following operational flow. The 
method for estimating a loop composition of a subscriber loop in terms of loop parameters 
Xi , X2 , . . . , Xi , . . . , Xn , the loop having a frequency-domain response 
H(o, Xi , X 2 , . . . , Xi , . . . , X N ) for the loop parameters, includes the following steps: 
Step 305 ~ determining a range for each loop parameter Xi , 
Step 3 10 - for each loop parameter Xi , generating a frequency-domain loop 
parameter function F Xi (#>) wherein 

F Xi (®W 1 ■» J »■ lX i H(^X 1 ,X 2? ... ? X i? ... 5 X N )dX 1 dX 2 ...dXi...dX N? 
x x x 2 x { X N 

Step 3 15 - generating a loop kernel k((D,fi) for all loop parameters wherein 

k(a>J) = J J ... J H(^X 1 ,X 2? ...,X N ) HGff,Xi,X 2 ,...,X N ) dX!dX 2 ...dX N , 

x x x 2 x N 

Step 320 — generating a parameter response function g { (fi) f° r eac h 1°°P 
parameter from the integral relation F Xi (&) = J k(oo, J3) g ; 

Step 325 ~ energizing the loop from a measurement end with an energy source, 
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Step 330 measuring a response signal H R (o?) = H(a> 3 Xi ) X2 5 ...,Xi,... > X N ) 

for the loop at the measurement end, and 

Step 335 - directly determining each loop parameter Xi from the integral relation 

Xi = iKMsMdP- 



Example in Terms of Specific Subscriber Loop Parameters 

The method for estimating a loop composition of a subscriber loop in terms 
of loop parameters such as the length of a cable section, denoted L, and the two possible 
gauges for the cable section, denoted Gl and G2, the loop having a frequency-domain 
response H(<2), L, G1,G2) for the loop parameters, includes the following steps: 

(a) determining a range for each loop parameter L, Gl, and G2, 

(b) for each loop parameter L, Gl and G2, generating a frequency-domain 
loop parameter function 

FlOW J J L H(<2, L, Gl , G2) dLdGl dG2 , 

L Gl G2 

Foi(*>) = M f Gl H(6>,L,G1,G2) dLdGldG2, and 

L Gl G2 

F G2 O) « J I i G2 H(a>, L, Gl, G2) dLdGldG2 

L Gl G2 

(c) generating a loop kernel k(a),fi) for all loop parameters wherein 
Ha> 9 fi) =! J J H(*>,L,G1,G2) H(#L,Gl,G2)dLdGldG2, 

L Gl G2 

(d) generating parameter response functions g L (/0, g G i(/?)> an ^ %gi(P) 
from the integral relations 
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F L (tf>) = Jk(»,£)g £ (040, 

ft 

Foi(©) = JktoflgoiGW.and 

P 

P 

(e) energizing the loop from a measurement end with an energy source, 

(f) measuring a response signal Hr(o) = H(<»,Lr,G1r,G2r) 

for the loop at the measurement end where Lr , G1r , and G2r are in or proximate to each 
range specified in step (a), and 

(g) directly determining each loop parameter L, Gl, and G2 from the 

integral relations 

L = fH R (/?)g L (/W, 
ft 

Gl= jH R (/?)g G1 (>W,and 
ft 

G2= jH R (^)g G2 (W. 

ft 

Although various embodiments which incorporate the teachings of the 
present invention have been shown and described in detail herein, those skilled in the art 
can readily devise many other varied embodiments that still incorporate these teachings. 
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What is claimed is: 



1 1. A method for determining configuration parameters describing a physical 

2 system, the method comprising the steps of 

3 measuring an output signal from the system in response to an input signal, 

4 the output signal being related to the configuration parameters by a linear operator, and 

5 directly reconstructing each of the configuration parameters by applying a 

6 prescribed mathematical algorithm to the output signal. 

1 2. The method as recited in claim 1 wherein said step of directly reconstructing 

2 includes the step of computing a configuration parameter function. 

1 3. The method as recited in claim 1 wherein said step of directly reconstructing 

2 includes the step of computing a configuration kernel 

1 4. The method as recited in claim 1 wherein said step of directly reconstructing 

2 includes the step of computing a configuration parameter response function for each of the 

3 configuration parameters. 

1 5. A method for estimating a loop composition in terms of loop parameters 

2 representative of the loop composition comprising the steps of 

3 energizing the loop from a measurement end with an energy source, 

4 measuring a response signal from the loop at the measurement end, 

5 wherein each of the loop parameters is related to response signal by a linear operator, and 



-24- 



Lechleider-Peoples Case 3-5 



1 directly reconstructing each of the loop parameters by executing a 

2 prescribed mathematical algorithm, determined with reference to the linear operator, on 

3 the response signal 

1 6. The method as recited in claim 5 wherein said step of directly reconstructing 

2 includes the step of computing a loop parameter function. 

1 7. The method as recited in claim 5 wherein said step of directly reconstructing 

2 includes the step of computing a loop kernel. 

1 8. The method as recited in claim 5 wherein said step of directly reconstructing 

2 includes the step of computing a parameter response function for each of the loop 

3 parameters. 

1 9. A method for estimating a loop composition of a subscriber loop in terms of 

2 loop parameters Xi , X 2 , . . . , Xi , . . . , X N > the loop having a frequency-domain response 

3 H(#>, Xi , X 2 , . . . , Xi , . . . , X N ) for the loop parameters, the method comprising the steps of 

4 (a) determining a range for each loop parameter Xi , 

5 (b) for each loop parameter Xi , generating a frequency-domain loop 

6 parameter function F Xi (#>) wherein 

7 Fxi(*>) = f f ... J ... IXiH(6>,Xi,X2,...,Xi,...,X N ) dXidX 2 ...dXi...dX N , 

XiX 2 x { X N 

8 (c) generating a loop kernel k(a>,/J) for all loop parameters wherein 
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1 k(&,/3) = J J ... J H(fi) J X 1 ,X 2? ... ? X N ) H(#Xi,X 2 ,...,X N ) dXidX 2 ...dX N , 

2 (d) generating a parameter response function g { (J3) for each loop 

3 parameter from the integral relation F Xi (<o) = J k(a>, p) g, (fl)dfi , 

P 

4 (e) energizing the loop from a measurement end with an energy source, 

5 (f) measuring a response signal H R (a>) = H(<5>, Xi , X 2 , . . . , Xi , . . . , Xn) 

6 for the loop at the measurement end, and 

7 (g) directly determining each loop parameter Xi from the integral relation 

8 Xi = iHRG^giOW. 

1 10. The method as recited in claim 9 wherein step (e) includes the step of 

2 computing the inverse of k(co,fl) . 

1 11. The method as recited in claim 9 wherein step (e) includes the step of 

2 computing the inverse of k(a>, P) using singular value decomposition. 

1 12. The method as recited in claim 1 1 wherein step (f) includes the step of filtering 

2 noise from the response signal. 

1 13 . A system for generating the loop composition in terms of loop parameters 

2 representative of the loop composition comprising 

3 a source of waves for energizing the loop from a measurement end, 
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1 a detector for detecting a response signal from the loop at the measurement 

2 end, wherein each of the loop parameters is related to response signal by an integral 

3 operator, and 

4 a reconstructor for directly reconstructing each of the loop parameters by 

5 executing a prescribed mathematical algorithm, determined with reference to the integral 

6 operator, on the response signal. 



1 14. The system as recited in claim 13 wherein said reconstructor includes a 

2 processor for computing a loop parameter Sanction. 

1 15. The system as recited in claim 13 wherein said reconstructor includes a 

2 processor for computing a loop kernel. 

1 16. The system as recited in claim 13 wherein said reconstructor includes a 

2 processor for computing a parameter response function for each of the loop parameters. 
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ABSTRACT OF THE DISCLOSURE 

A method for the direct estimation of configuration parameters defining the 
make-up of a physical system, such as cable lengths and gauges in the case of telephone 
subscriber loops, from boundary or end-point measurements of the physical structure, 
such as the impedance or return loss response of a subscriber loop, due to a probing 
source. 
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FIG. 2A 

^ 205 

DETERMINE LOOP CONFIGURATIONS THAT ARE EXPECTED 





^-210 


SELECT A COLLECTION OF THESE CONFIGURATIONS 


1 ^-215 


DETERMINE FREQENCY RESPONSES FOR THIS COLLECTION 




^ 220 

r r 


STORE ARRAYS REPRESENTING THESE RESPONSES 




S~ 225 


FORM AND STORE AN ARRAY OF ALL PAIRS 
OF FREQUENCIES FROM THE SET OF FREQUENCIES 
FOR WHICH MEASUREMENTS ARE MADE 




230 

r f 


FOR EACH PAIR OF FREQUE 
MULTIPLY THE RESPONSES FO 
STEP 210, SUM THE PRODUC 
AND STORE RES 


ncies stored in step 225, 
r the configurations from 
:ts over the collection, 
iult as a matirx 



r 



235 



INVERT THIS MATRIX AND INDEX THE INVERTED MATRIX 
BY PAIRS OF FREQUENCIES CALLED 
INPUT AND OUTPUT FREQUENCIES 



f GO TO^N 
V ^FIG. 2B J 



FROM 
FIG. 2A 



FOR EACH RESPONSE OF STEP 215, MULTIPLY THE REPONSE BY 
THE CONFIGURATION REPRESENTED BY THE RESPONSE, AND 
SUM THE PRODUCTS OVER ALL CONFIGURATIONS OF STEP 210 



I 



245 



MULTIPLY EACH ELEMENT OF THE MATRIX OF STEP 235 BY THE 
RESULT OF STEP 240 WITH THE SAME ELEMENT 
AS THE OUTPUT FREQUENCY 



I 



^ 250 



SUM THE RESULT OF STEP 235 OVER ALL OUTPUT 
FREQUENCIES FOR EACH INPUT FREQUENCY 
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255 



STORE THE RESULT OF STEP 250 AS A G-ARRAY 
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260 



MEASURE THE SYSTEM TO DETERMINE 
A MEASURED RESPONSE 



265 



FOR EACH MEASUREMENT FREQUENCY, MULTIPLY THE 
MEASUREMENT BY THE SUITABLE VECTOR IN THE G-ARRAY 



I 



270 



SUM THE RESULT OF STEP 265 OVER ALL MEASUREMENT 

FREQUENCIES 



FIG. 2B 
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DETERMINE A RANGE FOR EACH 
LOOP PARAMENTER 



310 



FOR EACH LOOP PARAMETER, GENERATE 
A LOOP PARAMETER FUNCTION 



I 
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GENERATE A LOOP KERNEL FOR ALL 
LOOP PARAMETERS 



I 



320 



GENERATE A PARAMETER RESPONSE 
FUNCTION FOR EACH LOOP PARAMETER 



r 



325 



ENERGIZE THE LOOP FROM A MEASUREMENT 
END WITH A SOURCE 



330 



MEASURE A RESPONSE SIGNAL FOR THE 
LOOP AT THE MEASUREMENT END 



335 



DIRECTLY DETERMINE EACH PARAMETER 
FROM AN INTEGRAL RELATION BETWEEN 
THE RESPONSE SIGNAL AND EACH 
PARAMETER REPONSE FUNCTION 



FIG. 3 
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IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 
Declaration and Power of Attorney 



As a below named inventor, I hereby declare that: 

My residence, post office address, and citizenship are as stated below next to my 

name. 

I believe I am the original, first and sole inventor (if only one name is listed 
below) or an original, first and joint inventor (if plural names are listed below) of the 
subject matter which is claimed and for which a patent is sought on the invention entitled 
DETERMINING THE COMPOSITION OF PHYSICAL STRUCTURES FROM 
BOUNDARY CONDITIONS, INCLUDING THE COMPOSITION OF 
SUBSCRIBER LOOPS, the specification of which is attached hereto. 

I hereby state that I have reviewed and understand the contents of the above 
identified specification, including the claims and any amendments filed therewith. 

I acknowledge the duty to disclose information which is material to patentability 
as defined in 37 CFR 1.56. 

I hereby claim foreign priority benefits under 35 U.S.C. 1 19(a)-(d) or 365(b) of 
any foreign applications) for patent or inventor's certificate, or 365(a) of an PCT 
international application which designated at least one country other than the United States 
of America, listed below and have also identified below any foreign application for patent 
or inventor's certificate, or of any PCT international application having a filing date before 
that of the application on which priority is claimed: 

None 

I hereby claim the benefit under 35 U.S.C. 1 19(e) of any United States provisional 
application(s) listed below: 

Serial No. 60/158,156 filed October 8, 1999 

I hereby claim the benefit under 35 U.S.C. 120 of any United States 
applications), or 365(c) of any PCT international application designating the United States 
of America, listed below and, insofar as the subject matter of each of the claims of this 
application is not disclosed in the prior United States or PCT international application in 
the manner provided by the first paragraph of 35 U.S.C. 1 12, 1 acknowledge the duty to 
disclose information which is material to patentability as defined in 37 CFR 1.56 which 
became available between the filing date of the prior application and the national or PCT 
international filing date of this application: 



None 



I hereby declare that all statements made herein of my own knowledge are true 
and that all statements made on information and belief are believed to be true; and further 
that these statements were made with the knowledge that willful false statements and the 
like so made are punishable by fine or imprisonment, or both, under 18 U.S.C. 1001 and 
that such willful false statements may jeopardize the validity of the application or any 
patent issued thereon. 

As a named inventor, I hereby appoint John T. Peoples (Reg. No, 28,250) to 
prosecute this application and to transact all business in the Patent and Trademark Office 
therewith. 

It is respectfully requested that all written communications from the Patent and 
Trademark Office in connection with this application be addressed to: 

John T. Peoples 
14 Blue Jay Court 
Warren, New Jersey 07059; 

any telephone calls should be made to (908) 580-9816. 



Full name of first joint inventor: Joseph William LecMeider 



Inventor's signature ^ o^y^^L k/^lUs^^ — ■ 
Date Sy> l^oa o 

Residence: 4 Harding Terrace (Mendham Township) 
Morristown, Morris County, New Jersey 

Citizenship: United States of America 

Post Office Address: (as above) 



Full name of first joint inventor: John Terrance Peoples 

Inventor's signature / • U'<**A* / ^ 
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Residence: 14 Blue Jay Court 

Warren, Somerset County, New Jersey 

Citizenship: United States of America 

Post Office Address: (as above) 
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